arXiv:1505.07195vl [nucl-th] 27 May 2015 


Dilepton production spectrum above Tc with a lattice quark propagator 
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The dilepton production rate from the deconfined medium is analyzed with the photon self¬ 
energies constructed from quark propagators obtained by lattice numerical simulation for two values 
of temperature T = l.STc and STc above the critical temperature Tc- The photon self-energy is 
calculated by the Schwinger-Dyson equation with the lattice quark propagtor and a vertex function 
determined so as to satisfy the Ward-Takahashi identity. The obtained dilepton production rate at 
zero momentum exhibits divergences reflecting van Hove singularity, and is significantly enhanced 
around w ~ T compared with the rate obtained by the perturbative analysis. 


I. INTRODUCTION 

Ultra-relativistic heavy ion collisions are the unique 
method to produce the deconhned medium experimen¬ 
tally on the Earth mm- Various observables are mea¬ 
sured in the experiments [1] to reveal properties of the de¬ 
confined medium and a variety of phenomena which come 
into play during the time evolution of the hot medium. 
Among these observables, the dilepton production yield 
has a characteristic feature that the yield provides us di¬ 
rect signals from the primordial deconfined medium [3], 
because dileptons once produced in the hot medium do 
not interact with and pass through the medium owing to 
their colorless nature. 

The dilepton yield observed experimentally consists of 
the sum of the dilepton production in each stage of the 
time evolution of the hot medium. The dilepton pro¬ 
duction in heavy ion collisions is roughly classihed into 
three processes except for final state hadronic decays. 
The first one is the hard process, in which dileptons are 
produced by scatterings of hard partons in the colliding 
nuclei. The second and third ones are thermal radiations 
from the deconfined medium and confined medium, re¬ 
spectively. The dileptons in low invariant mass region 
are usually expected to be dominated by these thermal 
radiations, while dileptons from the hard process have 
relatively high transverse momenta and large invariant 
masses. Experimental results on the dilepton production 
yield are usually compared with the baseline called cock¬ 
tail, which is the production yield estimated from the 
observed hadron abundances and their branching ratios 
into a dilepton pair. If there is no dilepton production 
from hot medium, dilepton production yield should be 
consistent with the cocktail result. 

At Relativistic Heavy Ion Collider (RHIC) at 
Brookhaven National Laboratory, e“'"e“ pair produc¬ 
tion yield is measured by two collaborations, STAR and 
PHENIX [H [5]. Both of these collaborations reported 
that the pair production yield measured in Au-l-Au colli¬ 
sions at y/sNN = 200 GeV has significant enhancements 
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in the low invariant mass (m) region compared with the 
cocktail HU], although the descrepancy on the magni¬ 
tude of the enhancement between the two experimen¬ 
tal groups has not been settled down yet. The result of 
PHENIX Collaboration shows that the yield at to ~ 500 
MeV is about one order larger than the cocktail |3]. 

When one performs an estimate of the dilepton pro¬ 
duction yield, one first calculates the dilepton produc¬ 
tion rate per unit time and unit volume from a static 
medium. The dilepton production yield is then given by 
the space-time integral of the rate from each volume el¬ 
ement of the medium. The dilepton production rate of 
a static medium is proportional to the imaginary part 
of virtual photon self-energy [2H9]. When the temper¬ 
ature (r) is asymptotically high, the photon self-energy 
can be calculated perturbatively. Using the hard thermal 
loop (HTL) resummed perturbation theory [101 E] > the 
dilepton production rate was calculated in Ref. [12] for 
lepton pairs with zero total three-momentum, and the 
result was extended in Ref. m to nonzero momentum. 
It is, however, nontrivial whether or not such pertur¬ 
bative analyses well describe the production rate from 
the deconfined medium near the critical temperature 
which has turned out to be a strongly-coupled system 
Moreover, it is known that the perturbative anal¬ 
yses in Refs. [121 IE] are modihed by proper inclusion 
of higher-order terms M- The analysis of higher order 
terms, however, is complicated and it is still under debate 
whether the scheme is valid for the whole kinetic region 
mum. Eor the description of the dilepton production 
rate in the deconfined phase near Tc, therefore, it is desir¬ 
able to evaluate the rate without resort to perturbative 
methods. In particular, the large enhancement observed 
at PHENIX [4| suggests a possibility that the dilepton 
production from the strongly-coupled medium above Tc 
has a large enhancement compared with the perturbative 
results used in the previous analyses |16j . 

There are several attempts of non-perturbative anal¬ 
yses for the dilepton production rate with lattice QCD 
miim- When one investigates the rate on the lattice one 
must take an analytic continuation from the imaginary- 
time correlator computable on the lattice to the real-time 
photon self-energy. This procedure, however, is an ill- 
posed problem, since the information on the imaginary- 
time correlators for discrete imaginary-time points ob- 




2 


tained on the lattice is insufficient to reconstruct the con¬ 
tinuous real-time function by itself [19]. In Ref. [I8|, an 
ansatz for the spectral function is introduced to avoid 
this problem. An alternative way is to use the Baysian 
analysis such as the maximum entropy method |19j . The 
lattice correlator, however, is insensitive to the structure 
of the spectrum in the low energy region |20) . The esti¬ 
mate of the low-energy spectrum on the lattice, therefore, 
is a difficult problem with the presently existing method¬ 
ologies. 


The exact non-perturbative photon self-energy can 
be calculated by the Schwinger-Dyson equation (SDE) 
if we have the full quark propagator and the photon- 
quark vertex function. Recently, an analysis of the non- 
perturbative quark propagator above Tc was performed 
on the lattice in the quenched approximation in Lan¬ 
dau gauge In this series of analyses, the quark 

propagator was analyzed with the two-pole ansatz for the 
real-time propagator in order to carry out the analytic 
continuation. This pole ansatz was employed motivated 
by the study of fermion propagator at nonzero tempera¬ 
ture [231 [33] . It has been shown [3TH33] that this simple 
ansatz can reproduce the quark correlator obtained on 
the lattice over a rather wide range of bare quark mass 
and momentum. It is therefore expected that the ob¬ 
tained quark propagator well grasps the gist of the non- 
perturbative nature of the quark propagator. 


The purpose of the present study is to analyze the 
dilepton production rate using this quark propagator. 
We construct the SDE with the quark propagators ob¬ 
tained on the lattice in Ref. [33] and with the vertex 
function constructed so as to satisfy the Ward-Takahashi 
identity (WTI). Our formalism, therefore, fulhlls the con¬ 
servation law of electric current. In this analysis, we show 
that the obtained dilepton production rate exhibits an 
enhancement of one order or more compared with the 
one from the free quark gas. Compared with the pertur¬ 
bative result in Ref. [T3], our result has a qualitatively 
similar behavior at low m region, while it exhibits an 
enhancement around m :^T owing to van Hove singular¬ 
ity. The effect of the vetex correction is also discussed in 
detail. 


The outline of this paper is as follows. In the next 
section we introduce the SDE for the photon self-energy 
and its components, lattice quark propagator and a vet- 
erx function satisfying the WTI. In Sec. m we then solve 
the SDE and present the form of the dilepton production 
rate in our formalism. The rate without the vertex cor¬ 
rection is also calculated in this section. We then present 
the numerical result in Sec. lYJ The hnal section is de¬ 
voted to a short summary. 



FIG. 1. Diagramatic representation of the Schwinger-Dyson 
equation for the photon self-energy, Eq. it- The shaded cir¬ 
cles represent the full quark propagator and the full vertex 
function. 


II. SCHWINGER-DYSON EQUATION FOR 
PHOTON SELF-ENERGY 


A. Schwinger-Dyson equation 


As dileptons are emitted from decays of virtural pho¬ 
tons, the dilepton production rate from a medium per 
unit time and unit volume is related to the retarded self¬ 
energy of the virtual photon as [^[S] 


dr 

duid^q 


a 1 1 

12^* 0/3w _ I 




( 1 ) 


at the leading order of hne structure constant a with 
= uP'—q^ and the inverse temperature /3 = I/T. With 
the SDE in Matsubara formalism, the exact photon self¬ 
energy is given by the full quark propagator S'(P) and 
the full photon-quark vertex T^{P + Q,P) as 


{iwm, q) = -J2 ^ [ 7 ^ 

TrcTro [S{P)^^S{P + Q)T,(P + Q, P )], (2) 


where uJm = 2'KTm and Vn = (2n -I- I)7rT with inte¬ 
gers m and n are the Matsubara frequencies for bosons 
and fermions, respectively, = {iVmP) is the four- 
momentum of quarks, and Cf is the electric charge of 
a quark with an index f representing the quark flavor. 
The color, flavor, and Dirac indices of S[P) are sup¬ 
pressed for notational simplicity. Trc and Trp denote 
the trace over color and Dirac indices, respectively. We 
note that since we take Landau gauge in this calcula¬ 
tion, off-diagonal elements in color space disappear. As 
a result, the trace over the color indices gives a factor 3 
in Eq. ([^. Equation ([^ is graphically shown in Fig. 
in which the shaded circles represent the full propagator 
and vertex function. The retarded photon self-energy is 
obtained by the analytic continuation. 

In the following, we consider the two-flavor system with 
degenerate u and d quarks, in which ef = 5e^/9. In 
this study we also limit our attention to the q = 0 case. 





B. Lattice Quark Propagator and Spectral 
Function 
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In the present study, we use the quark propagator ob¬ 
tained on the quenched lattice in Ref. [23] as the full 
quark propagator in Eq. ([^. In this subsection, af¬ 
ter a brief review on the general property of the quark 
propagator we describe how to implement the results in 
Ref. |23| in our analysis. 

On the lattice with a gauge fixing, one can measure 
the imaginary-time quark propagator, 

S^t^{T,p) = J (4) 

where x) is the quark field with the Dirac index /i. 
Here, r is the imaginary time restricted to the interval 
0 < T < /3. For the moment, the Dirac indices fi and v of 
the quark propagator are explicitly shown. The Fourier 
transform of the quark correlator, 

S^y{ii'n,p)= / (5) 

is written in the spectral representation as 

Sij,^{wn,p) = - [ dv' ( 6 ) 

7-00 

with the quark spectral function p^u{i'',p). The spec¬ 
tral function is related to the imaginary-time correlator 
Eq. Q as 

roo g(l/2-r//3)/3y 

S^,{t,p) = (7) 



FIG. 2. Open sympols show the momentum dependence 
of the parameters i^+{p), V-{p), and Z-{p)l{Z+{p) + Z-{p)) 
obtained on the lattice in Ref. [23]. The solid lines represent 
their interpolation obtained by the cubic spline method. The 
dashed line represents the light cone. 

On the lattice, one can measure the imaginary-time 
correlator Eq. for discrete imaginary times. To ob¬ 
tain the quark propagator one has to deduce the spec¬ 
tral function from this information. In Refs. pTM23| . the 
quark correlator in Landau gauge is analyzed on the lat¬ 
tice with the quenched approximation, and the quark 
spectral function is analyzed with the two-pole ansatz. 


In the deconfined phase in which the chiral symmetry 
is restored, the quark propagator anticommutes with 75. 
In this case, the spectral function can be decomposed 
with the projection operators A± (p) = (1 ± -fop ■ 7 ) /2 
as 


p{v,p) = p+(j/,p)A+(p)7o -bp-(z/,p)A_(p)7o, (8) 

with p = IpI, p = p/p, and 

P±{v,p) = ^Tro [p(i/,p)7oA± (p)]. (9) 

It is shown from the anticommutation relations of ^ and 
■0 that the decomposed spectral functions satisfy the sum 
rules. 


J diyp±{v,p) = 1. ( 10 ) 

Using charge conjugation symmetry, one can show that 
p±{i',p) satisfy [23] 


( 11 ) 


P-h(up) = Z+{p)6 {v - i/+(p)) -b Z_{p)d {v + v-ip )), 

( 12 ) 

where Z± (p) and v± (p) are the residues and positions of 
the poles, respectively. The four parameters, Z±{p) and 
(p) 1 are determined by fitting the correlators obtained 
on the lattice for each p. The two poles in Eq. ( [T^ at 
;/+(p) and v-{p), respectively, correspond to the normal 
and plasmino modes in the HTL approximation. In fact, 
the study of the momentum and bare quark mass, toq, 
dependences of the fitting paramters shows that 

the behavior of these parameters is consistent with this 
observation (211 [ 2 ^; for example, for large mg or p the 
residue of the plasmino mode (p) becomes small and 
the propagator approaches that of the free quark. The 
restoration of the chiral symmetry for massless quarks 
above is also checked explicitly on the lattice by mea¬ 
suring of the scalar term in the massless quark propagator 


In Fig. [2] w e show the fitting result of each parame¬ 
ter in Eq. (12) for massless quarks as a function of p for 
T = l.STc and 3Tc obtained in Ref. [23]. These analyses 
are performed on the lattice with the volume 128^ x 16, 
where both the lattice spacing and finite volume effects 


p±(ut) = Pt(-up)- 
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are found to be small [53]. In the upper panel, p depen¬ 
dences of v±{p), i.e. the dispersion relations of the nor¬ 
mal and plasmino modes, are shown by the open symbols. 
The vertical and horizontal axes are normalized by the 
thermal mass mx defined by the value of i^±{p) at p = 0: 
The value of mx obtained on the lattice after the extrap¬ 
olation to the infinite volume limit is tox/T = 0.768(11) 
and 0.725(14) at T = 1.5Tc and 3Tc, respectively [53] . 
The lower panel shows the relative weight of the plas¬ 
mino residue, Z-/ (Z+ -|- Z-). This figure shows that the 
weight becomes smaller as p increases, which indicates 
that the quark propagator for large p/T is dominated 
by the normal mode. A result similar to that shown in 
Fig-i is obtained by the Schwinger-Dyson approach for 
the quark propagator m- 

Although the lattice data are available only for dis¬ 
crete values of p, we must have the quark propagator 
as a continuous function of p to solve the SDE. For 
this purpose, we take the interpolation and extrapolation 
of the lattice data by the cubic spline method. From 
the charge conjugation symmetry one can show that 
diz_|_(p)/dp = —di/_(p)/dp, d^:/+(p)/dp^ = d^(p)/dp^, 
and Z+ (p) = (p) for p = 0 [1T| . These proper¬ 

ties are taken into account in our cubic spline interpola¬ 
tion. The lattice data are available only in the momen¬ 
tum range p/T < 4.7. To take extrapolations to higher 
momenta, we extrapolate the parameters using an expo¬ 
nentially damping form for Z- / (Z_|_ -|- Z -), 

Z_/{Z+ + Z_) = Re-^P, (13) 

and v±{p) are extrapolated by functions, 

r,^[p)=p + pfe-Ptp, (14) 

which exponentially approach the light cone for large p. 
The parameters R, a, and j3f' are determined in the cubic 
spline analysis. The p dependence of each parameter de¬ 
termined in this way is shown by the solid lines in Fig. 
We tested another extrapolation form by a polynomial, 
y±[p) = p+/p+/p^ + ■ ■ ■, but found that it hardly 
changes the dispersion relation. Finally, we fix 


Z+ + Z_ = l 


(15) 


throughout this paper to satisfy the sum rule Fq. (10 1 . 

With the two-pole form of the spectral function 
Fq. (j^, the quark propagator reads 


S{ivn,p) = 5 '+(zi/„,p)A+(p) 7 o -k S'_(zi/„,p)A_(p) 7 o, 

= X! Ss(,Wn,p)^s{pho, (16) 


S — ± 


propagator is given by 

S~^{Wn,P) = ^7^i^^n,Pho^s{p), (18) 

s—± 


with 


Ss ^{Wn,p) 


and 


{Wn - SV+{p)){Wn + SV-{p)) 

iu^-sE{p) ’ ^ ^ 


E[p) = -Z+{p)u_{p) + Z_{p)u+{jp). (20) 

Note that the inverse propagator has poles at = 
±A(p). These poles inevitably appear in the multipole 
ansatz, because the propagator Fq. 0 has one zero 
point in the range of u) surrounded by two poles. The 
form of the inverse propagator Fq. 0 will be used in 
the construction of the vetex function. We will see that 
the poles at = ±i?(p) give rise to additional terms in 
the dilepton production rate. 


C. Vertex Function 


The SDE, Fq. requires the full photon-quark ver¬ 
tex r^(P -|- Q, P) besides the full quark propagator. So 
far, the evaluation of r^(P -k Q^P) on the lattice at 
nonzero temperature has not been performed to the best 
of the authors’ knowledge. In the present study, we con¬ 
struct the vertex function from the lattice quark propa¬ 
gator respecting the Ward-Takahashi identity (WTI) as 
follows. 

The gauge invariance requires that the vertex function 
must fulfill the WTI 


EQ,P) = S-\P + Q) - S-\P), ( 21 ) 


with the inverse quark propagator S~^{P). For q = 0, 
the temporal component Tq is completely determined 
only by this constraint as follows. First, in this case 
q ■ r should vanish provided that T^ {i = 1,2,3) are not 
singular at q = 0. Then, by substituting q • F = 0 in 
Fq. (21) one obtains 


To{iojm + iun,p;ivn,p) 

= nh +iVn,P) - S~'^{wn,p)] ■ ( 22 ) 

On the other hand, the spatial components T^ cannot 
be determined only with Fq. (21) [55]. In the present 


study, we employ an approximation to neglect the q de¬ 
pendence of ro(fa;m+*^'ra,P+q; ^n, p) at q = 0 . In other 
words we assume that 


where 


Ss{wn,p) 


Z+{p) ^ Z_{p) 

iVn - SU+ (p) iVn P SV_ (p) ' 


(17) 


The symbols s = ± on the right-hand side are under¬ 
stood as the numbers ±1. Correspondingly, the inverse 


dTo{iu;m + iVn,p + q;iVn,p)/dqi\q=o = 0. (23) 


With this approximaiton and Eq. (21), one obtains 

^ Ej (iiOrfi ilZji ^ p Qi 1 p') 

= pivn,p + q) - S~'^{iu:m + iVn,p)- (24) 
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By taking the leading-order terms in q on the both sides, 
one has 

dS-^ 




-{iuJm + iVn,P) 


E dS^^{iuJm + iVn,p) . . , 
-HZi-70 A. (p) 


S=± 


dp^ 


S's + iVn,P)lO 


s=± 


dAJj^ 

’ 


(25) 


where in the second equality, we used Eq. (18). 


We note that there is no a priori justification of 
Eq. (23). By expanding Eg with respect to q at q = 0, 


one obtains 


iuj^To{iUJjn+Wn,P+ q\Wn,P) 

= + ivn,p) - S~^{wn,p) 

-fq • p-ioA{iuJm + iVn,P^) 

+ {q ■ p){p ■ l)B{iuj m H“ 5 ? P ) 

-bq • -yCiiuJm + ivn,iVn,P^) + 0{q^), 


(26) 


where A, B, and C are unknown functions. Our approx¬ 
imation corresponds to neglecting these functions. Al¬ 
though these functions do not affect Eq. (22) at q = 0, 
the corresponding terms appear in Eq. (25) when these 
functions are nonzero. The determination of the non- 
perturbative form of the photon-quark and gluon-quark 
vertices is generally difficult, and various approximations 
have been employed in studies of the SDE [5SHSI1- It 


should be emphasized that the vertex functions, Eqs. (22) 


and (25), satisfy the WTI and thus are advantageous in 


light of the gauge invariance among various ansatze on 
the vertex function. Eq, Eq. (22), is the same as that 
obtained in Ref. [28j since it is uniquely determined only 
from the WTL On the other hand, E^ differ from the 
ones in Ref. [28], even when only the longitudinal part 
in Ref. [2HI is concerned. Introduction of the functions 


given in Eq. (26) fills this difference. 


III. DILEPTON PRODUCTION RATE 

The goal of the present study is to obtain the dilepton 
production rate with the lattice quark propagators and 
the vertex function discussed in the previous section. In 
this section, however, before the analysis of the full ma¬ 
nipulation we first see the dilepton production rates in 
simpler cases: (1) free quark gas in Sec. Ill A and (2) 
case with the lattice quark propagators but with the bare 
vertex function in Sec. |III B| The full analysis is then pre¬ 
sented in Sec. Imcl 


A. Free quark gas 

The photon self-energy for the massless free quark gas 
is obtained by substituting the free quark propagator 
S{ivn,p) = l/(ir'n7o — P ■ l) and bare vertex function 
r,x = 7 ^ into Eq. (§. The result of the dilepton pro¬ 
duction rate for the massless two-flavor case is given by 


dV 

dcod^q 


q=0 


5q!^ 

367r'^ 


(/(l))^ 


(27) 


where /(w) = l/(e^‘^ -I- 1) is the Fermi distribution func¬ 
tion. 


B. Dilepton production rate without vertex 
correction 

Next, we calculate the photon self-energy and the 
dilepton production rate with the lattice quark propaga¬ 
tors Eq. (16) but with E^ = 7 ^. The photon self-energy 
obtained in this way, of course, does not fulfill the gauge 
invariance. The result obtained here, however, is helpful 
in understanding the effect of the vertex correction, i.e. 
the role of the WTI. 

The photon self-energy with the bare vertex is given 
by 


5e^ ^ 

n^i/(jwmj q') — 


d^p 


3 (27r)3 

Trials {WnJP)l^lS{iUJnl +iVn,p + q)lv\- 

(28) 


By substituting Eq. (16) into this formula, we obtain 

d^p 


0 ) = - 


{2i:f 

Ss{Wn,p)St{iuJm- 

s,t—zt 

X TrD[As7o7^At7o7,.]. 


ivn,p) 


The trace in Eq. (29) is calculated with 


TrD[As(p)At(p)] = 26 St, 
Tr£,[As(p) 7 o 7 iAt(p) 7 o 7 i] = 2ds-t + 2stpj, 
where it is understood that 

(5++ = 6— = 1, <5+_ = 5_+ = 0, 


Ss-t = (for t = ±)- 


(29) 

(30) 

(31) 

(32) 

(33) 


Substituting Eqs. (30), (31), and J2iPi = 1 irito Eq. (29), 
one obtains 


nC(*.;„,0)=^Ty] 


d^p 


3 (2^) 

Ss(iiyn,p)S-siiujjn + Wn,p), (34) 


s=± 
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where 

S-s{iVn + ii^m,p) = + iuJm,p) (for S = ±). 

(35), 


Using Eq. 0 and taking the Matsubara sum and the 
analytic continuation iojm to + ir), we obtain 


n;?-^(cc,0) 


—40a 
Stt , 


dp p^ 


Z+{p)Hl-2f{^+{p))) , Z.{p)Hl-2f{^.{p))) 


to — 2v+ (p) + irj 


to — 2v-{p) + ip 


+ 


2Z+{p)Z_{p) {f{v_{p)) - f+{v+{p))) \ 


00 — (p) + v_ (jp) + ip 




Taking the imaginary part of this result, we have 

A H/i /*oo 

Imn^’^(w, 0) =-^ dp p^ {Z+{pf5{oj - 2u+[p)) (1 - 2f{v+{p))) 

+ Z_{pf5{uj - 2v_{p)) (1 - 2f{v_{p))) 

+ 2 Z+{p)Z_{p)5{uj - v+{p) + v_{p)) {f{v_{p)) - f{v+{p)))} . 


(36) 


(37) 




FIG. 3. Photon production processes 


The imaginary part of the photon self-energy is the dif¬ 
ference between the annihilation and production rates of 
virtual photons in medium. The three terms in Eq. (37) 
represent different annihilation and production processes 
of a virtual photon. The first and second terms in 
Eq. (371 represent the production of a virtural pho¬ 
ton through the pair annihilation of two normal modes 
and two plasmino modes, respectively, which are dia- 


gramatically shown in Fig. (a), and their inverse pro¬ 
cesses. This can be checked from the arguments of the <5- 
functions and thermal factors in Eq. (371. The d-function 
in these terms represents the energy conservation during 
these processes, and the therrmal factor which is rewrit¬ 
ten as 


1 - 2f{oj) = (1 - f{oo)f - f{oo)\ (38) 

is the difference between the products of Pauli blocking 
effects and thermal distributions. The existence of the 
residues Z± (p) in these terms is similarly understood. 
The last term in Eq. ( [37| represents the Landau damp¬ 
ing between normal and plasmino modes, which is dia- 
gramatically shown in Fig. (b). Correspondingly, the 
thermal factor in this term can be rewritten as 


/(Wi) - f{uJ2) = /(Wi) (1 - /(W2)) - /(wa) (1 - /(wi)). 

(39) 


The ^-functions in Eq. (37) can be integrated out. The 
form of the dilepton production rate after the integration 
over p reads 




uj=2u+ (p) 


E 


pjZ-ipif 

2\dv^{pi)/dp\ 


(1 - 2f(v_{pi))) 


ui=2v-{pi) 


E 


‘^PiZ+{pi)Z_{pi) 
\d[v+{pi) - v-{pi)\/dp\ 


{f{u_{pi)) - f{v+{pi))) 


(40) 


UJ = V^{pi)-V_{pi) ^ 


where the momentum p in each term is given by the con¬ 
dition arising from the d-functions in Eq. (371. Each 


term can take nonzero values only when there exist mo¬ 
menta satisfying this condition for a given to. This gives 
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the condition for oj at which each term takes a nonzero 
value. For example, the first term takes nonzero values 
for u! > 2tot. Because the second and third terms can 
have multiple solutions of p for a fixed w, we represent 
this possibility by the sum over 1. It is also notable that 
each term in Eq. (401 is inversely proportional to the 


derivative of i^±ip) and i^+{p) — V-{p)', they come from 
the density of states of each modes. Accordingly, the 
dilepton production rate diverges when the derivatives 
vanish. Such divergence is known as van Hove singu¬ 
larity. In Sec. EYi we will see the appearance of such 


singularities in the dilepton spectrum. 

C. Dilepton production rate with vertex correction 

Now, let us calculate the dilepton production rate with 
the lattice quark propagtor Eq. (161 and the full vertex 
functions Eqs. ([22| and (25). 


When the full vertex function satisfying the WTI is 
used in Eq. ([^, the temporal component Hoo for q = 0 
vanishes. One can easily check this explicitly by substi¬ 
tuting E q (I 22 I into Eq. ([^. For_V^^^ H^^, by substitut¬ 
ing Eqs. (|16p and (25) into Eq. (^ one has 


^' 5g2 r ^3p 

, 0) — - T ^ ^ ^ S s{iT^njP')St{iT^n ’^^m-iP') 

i=l n s.t,u=± 


TA A A 1 , A » M 

2^ ---TrD[A„7o7*AtA„] -h TrD[A,,7o7iAt7o(7i - (P ' l)Pi)\ 


We substitute the following relations for the Dirac traces. 


Trp [As7o7jAtA„] = 2sp,6stStu, 

Trp [As7o7i^t7o7»] = 2^*-t + 2,stpj, 

Trp [As7o7iAt7o(p ■ 7 )] = ‘^Pi^st, 


and obtain 


^n,i(za;„,0) = f sSs{wn,p) 


fflnS's -I- iuJm,p) S-s{Wn + iUJm,p) 

dp p\ Ss{il^n+il^m,p) 


(41) 


(42) 

(43) 

(44) 


(45) 


where to obtain the first term we used 

a r .dS-\wn,p) d\nS-^{wn,p) 

^ - =-3^-. (46) 

Up to now, the calculation relies only on the decom¬ 
position Eq. (16), which is valid for the chiral symmetric 


quark propagator, and the form for the vertex Eq. (25). 
The result Eq. ( |45[ ) thus is valid for any form of the quark 
propagator S±(wn,p)- Now, we use the two-pole form of 
the quark propagator, Eqs. (|17| ) and (19). The term in¬ 
cluding d\nS~^/dp in Eq. (|45[) is then calculated to be 


a r- ^d\nS,^{Wn+iOJm,P) 

T sSs{lVn,p) g 

n ^ 


( 1 "l 

f dv+/dp 

dv-/dp 

dE/dp \ 

— SZ/+ Wn + SV- ) 

\iVr, + iuJm - SV+ 

iVn + iiOm + SV_ 

iVn + iuJm - sE) 


Z+dv-jdp 


Z-dv+/dp 


iUJm + Siy+ + SV- iuJm — SV+ — SV- 


{f{siy+) + f{si^-) - 1) 


Z+dE/dp Z_dE/dp 

U{sv+) - f{sE)) - — -—-— ifisiy-) - f{sE )), 


iuJm + 81^+ — sE 


iuJm — sv_ — sE 


( 47 ) 
































where the Matsubara sum over n is taken in the last equality. The remaining part of Eq. (451 is calculated as follows: 


T^sSs(iVn,p) ( 1 - 




Ss{wn + J P 


— 






Fi 


F, 


+ - 


F. 


iUn — i'L'n + SU- J \iVn + i^m — sE iUn + + i^m — SU- 

Zj^F\ ^ j, ^ , Z—F\ 


(48) 


iWm + SVj^ — sE 

Z+F 2 


ifisE) - fisv+)) 


f 2szy+ 
Z+F^, 


{f{-SV+) - /(Sl/+)) + T 


iUJm — sv_ — sE 

Z E2 


lUJra + Siy-f- — SV- 

Z-Fs 


{fisE)-f{-su.)) 


■ , -if{sv-) - fisi'+)) + ■ „ ifisv-) - fi-siy_)), 

lUJm, + s^+ — Siy- lUJm — 2,SV- 


(49) 


where 


Fi = - 
F2 = 

F3 = 


2E{v_ +E){v+-E) 
p{iy- — E){v+ + E) 
2v+{v+ - E){v+ - t/_) 
p{v++ E)(v++ v_) 
2v_{v^ — v_){v_ + E) 


p{i^++ - E) 

Here, each combination of i'± and E in the parantheses is set to become positive; this can be checked by the relation 
v+ > V- > —E > 0. 

Combining these results, Eq. (451 is calculated to be 
3 


^ ^ H^^ 0) — 


i=l 


lOe^ r <fp ^ f 2Zliy+fl 1 - 2/(i/+) 2Zliy_fi 1 - 2f(iy_) 

3 J (2Tr)^ \p{v++ E) iiUm+ ‘2si^+ p{iy-- E) iu)m +‘2sv- 

ClE - 2uj+uj f{v) - f{v+) 


+2z+z_n 


p{v+ + E){v- — E) iuJm — SV+ + SV- 

-[Z+—-Z ^ ~ ~ 

dp dp J iujjn + SZ/+ + sv- 

2Z+Z_E{v+ + v_ f dE\ / ^ f{E) - f{y+) ^ f{-E) - f{v_) 
p{v+ + E){v- — E) dp ) \ ^ iuJrn + — sE ioJm + SV- + sE 


(50) 


with H = i/+ — . 

By taking the analytic continuation iuJm OJ + irj, and taking the imaginary part, we obtain 


lmn«''^(cc,0)=-^ J dpp^ 


or) r y 2 y 2 

J + + (1 - 2f{v+)) S{oj - 2u+) + (1 - 2f{^_)) 6{oj - 212 .) 


p yv+ + E 


v_-E 


-z+z_ ~ -TT r (/(^"-) - f{^+)) -V+ + V-) 


{v+ + E){v- — E) 

(2Z+Z_E{v+ + v_f dE 
^ V p{'^+ + F){v_ -E) ^ dp 

X {Z+ (1 - f{-E) - f{v+))5{u: -V++ E) + Z_ {f{-E) - f{iy-)) S{uj - v_ - E)} 

+ (w —)■ —Uj). 


(51) 


Now let us inspect the physical meaning of each term in Eq. (511. From the i5-functions and thermal factors, one 
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finds that the two terms in the first line represent the pair 
creation and annihilation processes of normal and plas- 
mino modes, respectively. The second line corresponds 
to the Landau damping. These terms have correspond¬ 
ing counterparts in Eq. ( [3^ , although the coefficients of 
these terms are modified by the vertex correction. The 
term in the third line in Eq. (51) can be interpreted as the 


pair annihilation and creation of a normal mode and a 


plasmino one. This process does not appear in Eq. (37), 


and can manifest itself as a consequence of the vertex 
correction. We note that the similar process exists in the 
formula obtained in Ref. 1121 . In this way, the terms in 


the first three lines in Eq. (51) can be understood as the 


annihilation, creation, and scattering processes of quark 
quasi-particles. We also note that the Landau damping 
of two normal or two plasmino modes do not exist in 


Eq. (51), because such a process can exist only for w = 0 


at q = 0. 

On the other hand, one cannot give such interpreta¬ 
tions to the terms in the fourth and fifth lines in Eq. (51). 


From the (5-functions and the thermal factors, these terms 
seem to represent the decay and creation rates with a 
quasi-particle mode with energy ±if, which, however, 
does not exist in the quark propagtor Eq. (17). Mathe¬ 


matically, these terms come from the poles in the vertex 
function Eq. (25). The poles appear in the vertex func¬ 


tion via the WTI, Eq. (21) and the fact that the ana¬ 


lytic continuation of the propagator Ss{ii'n,p) gives zero 
points at energies ±E. As discussed in Sec. |IIB[ the zero 
in Sa{io,p) inevitably appears between the two poles in 
the two-pole form of the quark propagator Eq. 

Another remark on Eq. © is the sign of each term 
in Eq. (51). In Eq. (51), all terms are separately posi¬ 


tive definite for w > 0 except for the one in the third line, 
which becomes negative for sufficiently large lo. The neg¬ 
ative contribution of this term is, however, canceled out 
by the last term; we have checked that the sum of these 
terms is always positive. The total dilepton production 
rate for w > 0 therefore takes a positive value as it should 
be. 

In Fig. we show the dispersion relation of E{p) for 
T = 1.5Tc. For p = 0, E{p) vanishes because of chiral 
symmetry, while —E(p) approaches p for large p. The 
figure shows that —E{p) is a monotonically increasing 
function of p. The result for T = is qualitatively 
the same. In Fig. we also show the combinations of 
the dispersion relations appearing in the (5-functions in 
Eq. m, 


+ iy_, — V-, — E, v_ -{■ E. 


(52) 


From the figure, one sees that — E and u- -\- E are 
monotonically increasing and decreasing functions of p, 
respectively, starting from tot at p = 0. Also, — E > 
rriT and 0 < v- + E < tot are satisfied. These behaviors 
become transparent by rewriting these combinations as 

— E = Z+{v++ V-), (53) 

ly-+ E = (54) 


S’ 

q5 

c 

in 



0 0.5 1 1.5 2 2.5 3 3.5 4 


p/mj 


FIG. 4. Momentum dependences of various functions com¬ 
posed of iy±{p) and E. 


where we used Eq. (20). 


We finally comment on the limiting behaviors of 
Eq. ( [ 51 ] ). First, in our two-pole ansatz the quark propa¬ 
gator for massless free quarks is obtained by setting 


Z+{p) = 1, Z_{p) = 0, v+{p) = p. 


(55) 


Equation (51) thus should reproduce the photon self¬ 


energy of the free quark gas, when Eq. (55) is substituted. 


This can be explicitly checked as follows. By substituting 
Z- = 0, all terms including Z_ vanishes. Since E = —v- 
for Z- = 0, the third and fourth lines in Eq. (51) cancel 
out with each other without constraints on i/- (p). Only 
the first term in Eq. ( |51[ ) thus survives, which gives the 
free quark result. Second, our result on the dilepton pro¬ 
duction rate approaches the free quark one in the large lu 
limit, because the lattice quark propagator used in this 
study reproduces Eq. (55) at large momentum. This be¬ 


havior will be explicitly checked in the next section. 


IV. NUMERICAL RESULTS 


Now let us see the numerical results on the dilepton 
production rate obtained in the previous section. In 
Fig.[5l we present the ui dependence of the dilepton pro¬ 
duction rate for T — 1.5Tc. In the figure, we also plot the 
result without vertex correction, Eq. (37), together with 
the rates obtained by the HTL calculation m and the 
free quark gas, Eq. (27). The value of the thermal mass 
mx is taken from the the result obtained on the lattice 


m- 

Figure shows that the production rate with the lat¬ 
tice quark propagators has divergences at two energies, 
cu/mx = wi/mx — 1.1 and w/mx = a; 2 /m-x — 1-8. For 
uj < LOi, our result, as a whole, behaves similarly to the 
HTL one [H], i.e. it increases as u decreases, although 
our production rate is smaller than the perturbative one 
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colrrij 


FIG. 5. Dilepton production rate at zero momentum for T = 
l.STc. The result without vertex correction is also plotted. 
Thin lines represent the HTL and free quark results. 


for small lo. Near coi, however, it shows a prominent en¬ 
hancement and exceeds the latter. The region, where the 
large production rate is obtained, is located around rriT. 
Therefore, it can be possible that the production yield ob¬ 
tained by integrating this rate has the large enhancement 
below several hundred MeV, where the enhancement in 
the experimentally-observed dilepton spectra at RHIC 
Hi] exists. The rate has a discontinuity at ui = oji, 
and is significantly suppressed compared with Eq. (27) 
for wi < uj < UJ 2 - The rate has another discontinuity at 
to = UJ 2 , above which the rate is close to the free quark 
one. In the dilepton rate without vertex correction, one 
also finds two divergences at w = and UJ 2 , while the 
rate vanishes for wi < to <ui 2 - 



FIG. 6. Decomposition of the dilepton production rate. 


rately in Fig. In the figure, the rate coming from the 
pair annihilation of two normal modes (NN), two plas- 
mino modes (PP), and a normal and a plasmino modes 
(NP) are separately shown, together with those of the 
Landau damping between quasi-particles (LD) and pro¬ 
cesses including an E mode with a normal (NE) and a 
plasmino (PE) modes. From the figure, one finds that 
the two divergences at w = wi and UJ 2 come from the LD 
and PP rates, respectively. As discussed in Sec. [nT¥ 


these divergences come from van Hove singularity. The 
photon self-energy Eq. (51) is inversely proportional to 


derivatives of the dispersion relations, di'-{p)/dp and 
dlv+i^p) — iy-{p)}/dp. As shown in Figs. [^and|^ each 
of z/_ (p) and (p) — v- (p) has an extremum at nonzero 
p. Their values at the extrema are I'+ip) — v_{p) = wi 
and 2v-{p) = u) 2 - At these points the derivatives vanish. 
This leads to the divergences in the photon self-energy, 
and accordingly the dilepton production rate. 

Figure]^ also shows that each indivisual process is non¬ 
vanishing in a limited range of w. The range can be read 
off from the corresponding functions plotted in Fig. 
The NN and NP rates are nonvanishing for w > 2mT. 
On the other hand, the lower threshold of the PP rate, 
UJ = l>J 2 , is slightly lower than 2tot, because v-{p) has a 
minimum smaller than mx at nonzero momentum. The 
range of the Landau damping is also kinematically con¬ 
strained to a; < Wi. The NE and PE rates are nonvan¬ 
ishing for w > Tot and uj < tot, respectively. The NE 
rate gives rise to a nonzero value for uji < uj < UJ 2 - To 
the dilepton pruduction rate without the vertex correc¬ 
tion, only the NN, PP, and LD contribute and the rate 
vanishes for wi < w < 012 . 

For large w, the rate is dominated by the NN. This 
is a consequence of the fact that the quark propagator 
approaces the free quark one as p becomes larger. A 
glance at Fig. might give an impression that the NE 
rate also survives for large uj. Although not shown in 
Fig. [6| however, the NP rate takes a negative value for 
UJ > 2.4mT, and this term almost cancels out with the 
NE rate; see the discussion in Sec. |III C[ 

Next, let us address the behavior of the production 
rate in the w —^ 0 limit. The photon self-energy is identi¬ 
cal with the electromagnetic current-current correlation 
function, 

4[uj,p)= f 

(56) 

at the leading order in a. The low energy behavior of 
J^{uj,p) is related to electric conductivity a through the 
Kubo formula, 


1 1 ^ 

= - lim — J^{uj,0). 

6 UJ ^ * 


(57) 


i=l 


In order to understand these results in more detail, 
we show the contribution of each term in Eq. (51) sepa- 


Our result shows that X]i=i 0) approaches zero 

faster than uj in the w —> 0 limit, and thus the electric 
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conductivity vanishes. Incorporation of the width of the 
quasi particle modes, which is not included in the form 
of the quark propagator used in this study, may lead to 
nonzero a. 



ca/rrij 


FIG. 7. Dilepton production rate from STc deconfined phase. 

In Fig. we show the dilepton production rate at 
T = 3Tc. One sees that the behavior is qualitatively 
the same as the result for T = l.STc. In particular, there 
exists an enhancement of the rate around u ~ 1.5tot ow¬ 
ing to van Hove singularity. This indicates that the large 
enhancement of the dilepton production rate around 
w ~ Tot is a general result irrespective of T. The gap 
between wi and 0 J 2 is narrower than that at T = l.STc, 
because of the change of the dispersion relations v± (p) 
obtained on the lattice. One also finds that the rate takes 
a finite value at a; = 0, while it diverges for T = l.STc. 


This limiting behavior may depend on the way of the 
extrapolation of v± (p) to large momentum. 

V. SUMMARY 

In this study we have investigated the dilepton pro¬ 
duction rate using a quark propagator obtained on the 
lattice with a pole ansatz. The Schwinger-Dyson equa¬ 
tion for the photon self-energy is solved with the lat¬ 
tice quark propagator and the photon-quark vertex sat¬ 
isfying the Ward-Takahashi identity. The effect of the 
vertex correction is discussed by comparing the result 
with the calculation without vertex correction. Our nu¬ 
merical result shows that the dilepton production rate 
with the lattice quark propagators is larger by about one 
order or more compared with that with the free quark 
ones in the low invariant mass region. Compared with 
the HTL result, there exists a significant enhancement 
around uj ~ mx — T owing to van Hove singularity. 
This result is interesting since such a large enhancement 
in the deconfined medium near Tc may explain the ex¬ 
cess of the dilepton yield observed at PHENIX [3] in the 
low invariant mass region. To understand the effect of 
the enhancement of the dilepton rate to the experimen¬ 
tal result quantatively, the analysis with dynamical mod¬ 
els describing the space-time evolution of the hot medium 
and integration of the dilepton production rate is needed, 
which would be performed elsewhere. The comparison 
with other non-perturbative approaches such as Ref. |32] 
will be an interesting future work. 
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